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We study the transport in ultrathin disordered film near the quantum critical point induced by the 
Zeeman field. We calculate corrections to the normal state conductivity due to quantum pairing fluc- 
tuations. The fluctuation-induced transport is mediated by virtual rather than real quasi-particles. 
We find that at zero temperature, where the corrections come from purely quantum fluctuations, 
the Aslamazov-Larkin paraconductivity term, the Maki-Thompson interference contribution and 
the density of states effects are all of the same order. The total correction leads to the negative 
magnetoresistance. This result is in qualitative agreement with the recent transport observations in 
the parallel magnetic field of the homogeneously disordered amorphous films and superconducting 
two-dimensional electron gas realized at the oxide interfaces. 
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Introduction.- According to the microscopic BCS- 
theory [l[ magnetic field extinguishes superconductivity. 
In the absence of spin-orbit interaction there are two ba- 
sic mechanisms. The first one is diamagnetic effect asso- 
ciated with the action of the magnetic field on the orbital 
motion of electrons forming a Cooper pair. The second, 
paramagnetic mechanism, is due to Zeeman splitting of 
the states with the same spatial wave function but op- 
posite spin directions. In the former case, the estimate 
for the upper critical field follows from the condition 
H C 2^ 2 ~ $o, where $o = hc/2e is the flux quantum, 
£ = y/hDjA is the coherence length for the disordered 
superconductor with A being energy gap, and D diffusion 
coefficient. In contrast, Zeeman splitting destroys super- 
conductivity at the other critical field that follows from 
the condition gLl^BH z ~ A, where \ib = eh/2mc is the 
Bohr magneton and is renormalized giro-factor. The 
ratio between two fields is H Z /H C 2 ~ kpi 3> 1, where 
hp is Fermi momentum and £ is the elastic scattering 
length. Thus, in bulk systems, the suppression of su- 
perconductivity is typically governed by the first - dia- 
magnetic mechanism. The situation changes in the case 
of restricted dimensionality. For example, in the case 
of thin-film superconductor the above ratio changes to 
H Z /H C 2 ~ (kp£)(d/£), which can be small provided that 
film is thin enough d <C £/kp£, such that spin effects 
dominate. 

The scenario of paramagnctically limited superconduc- 
tivity has long history that goes back to pioneering works 
of Clogston and Chandrasekhar 0] . The first order phase 
transition from superconductor to paramagnct was found 
at the critical field approaching E z = \/2A at low tem- 
peratures. In practice, the measured film resistance fol- 
lows a hysteresis loop instead of a sharp first order 
transition. At low field E z < \/2A and zero tempera- 
ture the system is superconducting. With increasing the 
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FIG. 1: [Color online] Above the tricritical point T* the sec- 
ond order paramagnet to superconductor transition occurs 
along the (black) solid line obtained from Eq. {TJ). At T < T* 
this line becomes a supercooling part of the hysteresis, and 
the (blue) dashed line is its superheating part. The latter 
is obtained following Ref. [3] . The grey shaded area with the 
critical point (0, Ao) as its lowest corner bounded by the black 
dashed line marks the region of quantum fluctuations. 

field y/2A < E z < 2A the film is trapped in a super- 
conducting metastablc state. At fields exceeding the su- 
perheating threshold E z > 2A the film becomes normal. 
When the field is reduced back to zero, the normal state 
is metastable in the interval A < E z < \/2A In this 
paper we study the transport properties at the onset of 
transition to the superconductivity near the supercooling 
field Ef = A [see Fig. [J. This field E Z C (T) corresponds 
to the zero binding energy of a Cooper pair and can be 
determined from the standard equation [8| 

ln(T c /T c0 ) = ^(1/2) - Re ^(1/2 + i£f /4ttT c ) (1) 

similar to that in the theory of paramagnetic impuri- 
ties 0]. Here ip is the digamma function and T c q = 
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T C (H = 0) is the critical temperature in the absence of a 
magnetic field. The zero temperature solution of Eq. (JXJ) , 
E* C (Q) = A, defines the quantum critical point (QCP), 
which is premier interest of our study. 

Motivation- The renewed interest in the physics of 
paramagnetically limited superconductors is motivated 
by the rapid growth of its experimental realizations. Re- 
cent parallel magnetic field studies of two-dimensional su- 
perconducting systems were extended to much lower tem- 
peratures thus making it feasible to approach the limit 
of QCP. Tunneling spectroscopy of ultrathin Al and Be 
films revealed field-induced spin mixing and anomalous 
resonances in the density of states 0, EdJllL The lat " 
ter was successfully explained in theory 12l. 13 1. which 



emphasized the crucial role of superconducting pairing 
correlations in the paramagnetic state even far from the 
transition region. A surprising enhancement of supercon- 
ductivity by a parallel magnetic field, deduced from the 
transport measurements, was observed in ultrathin, ho- 
mogeneously disordered amorphous Pb films and the two- 
dimensional electron gas realized at the interface of oxide 
insulators LaA103 and SrTiOs In addition, pro- 

nounced negative magnetoresistance (NMR), concomi- 
tant with the enhanced T c , was reported. Although we 
do not dwell onto the issue of T r enhancement in these 



systems (see Ref. |15| for the recent theoretical propos- 
als), we show that transport anomalies, such as NMR, 
can be successfully addressed within BCS theory. 

The issue of NMR in superconductors, either near 
the QCP or near the parallel-field-tuned superconductor- 
insulator transition, was previously discussed in the lit- 
erature experimentally [it 



cally 



171 ] and attributed thcorcti- 



2G[ to the proliferation of superconductive fluc- 
21| . These studies emphasized mainly the or- 
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bital effect of a magnetic field on the pre-formed Cooper 
pairs. In this work we develop transport theory of param- 
agnetically limited ultrathin superconductors focusing on 
the quantum regime of zero temperature near the criti- 
cal Zeeman field. The regime of classical fluctuations was 
partially discussed in the early papers (22- 25|. 

Theory- In the vicinity of the transition transport 
properties of superconductors are governed by the fluc- 
tuation effects. These are famous paraconductivity phe- 
nomena introduced by Aslamazov and Larkin (AL) (26f . 
Maki and Thompson (MT) j27|, and also related den- 
sity of states (DOS) effects discussed first by Abrahams 
et al. [Hj]. We follow these classical papers and ap- 
proach the problem based on the diagrammatic pertur- 
bation theory. Note that the technique based on the 
time-dependent Ginzburg-Landau formalism applied for 
studying transport near QCP [29|, Ho| accounts correctly 
only for the classical part of AL-type contribution to the 
conductivity, but it misses completely the quantum zero- 
temperature corrections. Microscopic approach takes 
care of all the contributions including DOS part, resulting 
from the depletion of the normal state density of states 



by superconducting fluctuations, and also MT interfer- 
ence term (l8l - [20| . In fact, at T = where the cor- 
rections come from purely quantum fluctuations, these 
effects turn out to be of the dominant nature. In calcu- 
lations we assume diffusive limit, 
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Conditions ([2]) arc satisfied in many experiments 0, flo| . 

Within Kubo linear response formalism conductivity 
is obtained from a = —K R (oS)Jioj by analytic con- 
tinuation of the Matsubara current correlation kernel 
K{uj n ) = - J 1/T dTe tu " T (T T J(T)J(0)). This kernel can 
be conveniently presented as a sum of three contributions 
K = Kal + Kmt + Kdos- The general expression for 
the AL term reads (hereafter h — ks = 1): 

K AL {ui n ) = -e 2 T ^2 B Q,n k ^ n L Qflk L Qfl k +u> nl (3) 
where flk = 2irkT. The triangular vertex function 
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consists of two Cooperons 
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t(DQ 2 + \e„ - e m \ - iaE z sgn(e n - e„ 



(6) 

and an integral over the block of three Green's functions 
with G Pen = (ie„ — £p + (7i? 2 /2+sgn(£, 1 )/2T)~ 1 . Here we 
used notations: e m = 2irT(m + 1/2), £p = P 2 /2m — Ep, 
v P = dp£p, 0(e)-step function and sgn(e , )-sign function. 
Finally, propagator of fluctuating Cooper pairs in Eq. ([3]) 
is given by 
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When calcu- 
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lating B-vertex one should follow few basic steps [21| . 
i) To the leading order in the momentum transferred 



Q one can approximate G_ P 



+Q,-e„+Q k 



G 



P,-e n +Q k 



(vp ■ Q)(Gp_ s +Sl ) . ii) Furthermore, one can neglect 
Zeeman energy as compared to the inverse scattering 
time in the Green's functions [provided the condition of 
Eq. ((2J] and then completes P-integration in a standard 
way — > v J d£p J ^ £ -- in) Next is the fermionic 
Matsubara £ m -sum in Eq. Q , which can be found in the 
closed form with the result 
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iv) The remaining step of calculation is bosonic fife- 
sum followed by an analytical continuation iuj n — > lo. 
The latter are accomplished via the contour integration 
over the circle with two-brunch cuts at Imfi = 0, — lu u 
where the product of propagators in Eq. ^ has breaks 
of analyticity. After w-expansion of K^ l (uj) to the lin- 
ear order one finds for the AL conductivity correction 



a AL = a AL + a AL + a AL^ where 
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x{<9 w (Bqj 2iW ) 2 (Lqj 2 ) 2 - 9 w (Bqj 2 _ u ^y{Lq^y 

+d„[(B%?n- u ,J 2 " (^o,J 2 ]l^,ol 2 }- (11) 

The superscripts R/A in the vertex function and prop- 
agators stand for the retarded/advanced components 
while subscripts cl/q refer to classical/quantum. This 
convention comes form the observation that as T — > 
classical contribution vanishes while quantum remains fi- 
nite. 

We turn now to the derivation of the MT contribution 
whose response kernel is given by 



e 2 T L Q,nJ^Qji k 
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where 
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Jmt = ^2 v P v Q-P G P,s m+n G Q-p.n k _ 



X G P,e m G Q-P-e m +n k - (I 4 ) 

Momentum integration in the block of Green functions 
Jmt is done under the same approximations as in the 
case of AL term described above. According to the stan- 
dard convention 2l| we split now MT term into the so- 
called regular and anomalous contributions: 
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After the analytical continuation these translate into the 
ty cor: 

e 2 vD 



conductivity correction a"" " = cr^'eg ^ 
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In order to regularize logarithmically divergent momen- 
tum integral in the case of anomalous contribution we 
have introduced pair-breaking cutoff parameter T^. 

We finally discuss the density of states contribution to 
the conductivity. The latter is given by the similar to 
Eq. (fT2]) expression with 
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where 



S Q.?2fc, w „ = 2T E ( X Q,e m ,n k -e m ) 2j DOS, (18) 

J dos = Y j v 2 p(Gp em yGp em+Uri [G Q a _p nk _ em 
p 

I^E( G keJ 2 G Q V,o fc -eJ- (19) 
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After standard steps outlined above one arrives at the 



conductivity correction a DOS 
form 
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The equality between the two contribution in Eq. (|20b[) 
has parallels with the original fluctuation transport con- 
siderations at T — T c <C T, In the original near-T c prob- 
lem, the typical energy of diffusing pairs DQ 2 ~ T — T c is 
smaller than the thermal energy of quasiparticle ~ T. In 
our case, E z adds to the energy of pairs making it bigger 
than T . Correspondingly, unlike the near-T c case, the 
off shell energy of a pair, 2s ~ T. falls below the pair 
excitation energy set by E z . This causes a sign inversion 
of the energy denominator associated with the unbound 
intermediate state and the correction (|20bp turns to be 
positive. In general, derived above conductivity correc- 
tions are applicable at any field H and temperature T 
above the transition. In the following we discuss limiting 
case of interest. 

Results - It is convenient to regroup all contributions 
and present total conductivity correction as the a sum 
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of zero-tempcrature (So~ q ) and finite-temperature (5<tt) 
terms, namely 



6a(H, T) = 5cr q {H) + Sa T (H, T). 



(21) 



The first term here is determined by the quantum AL 
[Eqs. ([TP ]) -([n j) ] and DOS [Eq. (|20b|) ] contributions, and 
also regular part of the MT conductivity [Eq. (|16a|) ] . The 
remaining terms define Sot- The magnitude of Sa q de- 
creases monotonically with increasing field; this leads to 
a negative magnetoresistance at zero temperature. At 
finite temperature, based on how the quantum critical 
point is approached, there are several regimes that show 
different T and H dependencies, which should be ex- 
perimentally accessible. Below we focus on QCP only 
and extract the leading singularity in Sa q as the function 
of Zeeman field. Thermal contribution <5<tt and various 
crossover regimes will be discussed elsewhere (3l| . 

At zero temperature ±in -> ln[(DQ 2 ± iQ + 
io~E z ) I '4ttT] and the pair-propagator can be taken in the 
leading pole approximation 
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which is obtained from Eq. ([7]) under the conditions 
DQ 2 < A and \E C ± 0| < A . Here A = nT c0 /2 lE 
where InjE ss 0.57 is the Euler constant, and E c = 
\J ' E 2 — A§ . The branch cut of the propagator (due to the 
logarithmic structure) also contributes to 5a q but gives 
the sub-leading singularity. Within the same accuracy 
we compute vertex functions: 



m 



DQ 2 (DQ 2 ± itt)(DQ 2 ± ifi - 2iu), 

(23) 



0B£U 2 = %^{DQ 2 f (DQ 2 - 2iu). (24) 

All together this leads to the conductivity correction near 
the Zccman field-induced quantum critical point 
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which is obtained within the logarithmic accuracy. Equa- 
tion (f25|) is the main result of the paper. 

Discussions - The conceptual difference of our anal- 
ysis from the problem of fluctuation-induced transport 
close to T c is that unpaired particles, have finite exci- 
tation energy E z , see Eq. ([5]). As a result, the activa- 
tion probability of such pairs is suppressed exponentially 
cx exp(— E z /T). We argue that while in the standard 
case the real gapless pairs are only important in our case 
such pairs are always virtual. 

Let us illustrate this point taking AL correction as 
an example. Consider first standard case near-T c . In 



Eq. ([3]) the triangular vertex Eq. ([8]) can be estimated as 
Bq.xi,^ oc DQ x dU Q ^/d^l. Here rfg.n = L Q \ 2 + g~ Y is a 
particle-particle polarization operator with momentum Q 
entering in a DQ 2 — ifl combination. At small momenta 
we can take IIo ; f2 in the clean system. The imaginary 
part of the polarization operator Imll sa f dt;[n(— £ p + 
O)n(^) - h(-C p + n)h(Z p )]5{n - 20 = u(Cl/2) tanh 
where the particle and hole occupation numbers are 
n(e) = (1 + e e / T )~ x , n(e) = 1 — n(s). The real part, 
due to virtual pairs RcII w log | (D, 2 — T 2 )/u>^ | , is a famil- 
iar Cooper logarithm. The imaginary part contribution 
Bq.q.uj oc DQ x /T. In contrast, the real part contribu- 
tion vanishes at O = due to the particle-hole symmetry, 
v(Sl) = v. The expansion in ~ T — T c <£i T yields a 
correction small in the parameter (T — T c )/T c 1. 

In the presence of Zeeman field the situation is very 
different. The pair activation rate, Imll « v(yi)[n(u) /2 — 
E z /2) — n(uj/2 + E z /2)], gives exponentially suppressed 
contribution oc DQ x exp(—E z /T)/T. The real part, 
due to virtual pair excitation, can be obtained by the 
Kramers-Kronig relation, Rell « log |(f2 2 — E z )/uj^\. Its 
contribution to Bqsi.u is suppressed only algebraically 
oc DQ X T / E 2 . Unlike the standard case the virtual quasi- 
particles make a dominant contribution to the triangular 
vertex excitations. The algebraic suppression of vertexes 
is most pronounced in the case of the AL and is mani- 
fested in additional factors of DQ 2 , tt in Eq. ([2l?l) -([2"4" l) . 
which makes it logarithmic in E z /E c . Note that in the 
case of near-i? C 2 problem [18( the AL contribution is also 
suppressed due to the current matrix elements connect- 
ing adjacent Landau levels 

The regular MT and DOS contributions are propor- 
tional to a second derivative of the real part of the po- 
larization operator RcIIq.q. Since the latter is finite at 
fl = 0, these contributions are as singular as AL terms. 

We have checked explicitly that other contributions 
such as diffusion coefficient rcnormalization as well as 
contribution with only one or no Cooperon vertexes are 
either small or non-singular. Since the temperature can 
be set to zero in integrations over fast fermion degrees 
of freedom, the additional factors of r results in small 
prefactors tE z , tDQ 2 or rf2. 

Outlook- The spin-orbit scattering and finite thick- 
ness effects modify the fluctuation transport, due to the 
finite spectral weight in the particle-particle channel at 
zero frequency. Addition of a finite spin-orbit scatter- 
ing introduces a finite life time T -1 to the Cooperon. 
At lowest temperatures the superconductivity survives if 
this scattering is not too strong, L <g; E z with somewhat 
lower critical field. While E z approaches the supercooling 
transition from above the results obtained in the present 
paper are expected to cross over to a different regime at 
L a* E c . The finite film thickness affects the crossover 
in a similar way. All these relevant perturbations as well 
as the regime of close proximity to the supercooling line 
will be studied elsewhere 311. 
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